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We study the "mixed spin" isotropic ladder system having 5 = 1 spins on one leg and S = ^ 
spins on the other, with general-type exchange interactions between spins on neighboring rungs. A 
set of model Hamiltonians with exact ground states in the form of a certain matrix product wave 
function is obtained. We show that sufficiently strong frustration can lead to exotic singlet ground 
states with infinite (exponential) degeneracy. We also list a couple of rather simple models with 
nontrivial ground states, including a model with only bilinear exchange. 
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I. INTRODUCTION 

During the last few years, "mixed" one-dimensional 
(Id) quantum spin systems composed of two or more dif- 
ferent kinds of spin have drawn certain interest. Sev- 
eral Bethe-ansatz solvable models with singlet ground 
states have been foundplOexperimentally relevant mod- 
els (see, e.gp. jEjef. ^ of ferrimagnetic quantum chains 
were studiedutll via various numerica.1 a.nd analytical ait. 
proaches; recently, mixed spin chaing|13llj and laddersEil 
with singlet ground states were analyzed both numeri- 
cally and by means of the nonlinear sigma model tech- 
nique, and other interesting mixed-spin models prnth rich 
phase diagram were proposed and investigated.li3 

On the other hand, there has been a considerable 
progress in studying Id spin systems by means the 
so-called matrix product (MP) states technique MtIj MP 
states have proved to be a particularly useful tool for con- 
structing new models with exact ground states: ground 
states of the MP type were rfound for a large famibi. pf 
spin-1 and spin-| chainsEJtiJ and spin-i ladders Eja 
The simplest example of the MP state is the spin-1 va- 
lence bond stsXe (VBS) which is the ground state of the 
AKLT modelcil and is widely accepted as a convenient 
image of the Haldane-phase state. MP states were also 
successfully used for the variational study of the, Kpund 
state and,|elementary excitations of spin chaingHa^EZl and 
ladders.E3o However, the states which can be accessed 
via the MP approach have finite, typically rather short, 
correlation length. Therefore they are usually gapped 
(except for the case of ferr omagnpti r- ty p e situation with 
spontaneously broken symmetryllil'tj when the system 
has long-range order and is gapless due to the Goldstone 
modes) and cannot be applied for a description of the 
system behavior at any critical point. 

In this paper, we study the isotropic mixed-spin system 
which may be viewed as a ladder composed of 5 = 1 and 
S — ^ legs, with general bilinear, biquadratic, and six- 



spin interactions between neighboring rungs, as shown 
in Fig. 1^ The model is described by the Hamiltonian 
H = hn,n+i, where h is chosen in the following form: 

+ i J.(5lTi) + ^J'r{S2T2) + Jrf(xiS2) + ^^(^1X2) 

+ Ks{SiS2)^ + Ksr{S^S2){T^T2) 
+ Krr{SlTi){S2T2) + K m{t lS2){S 1T2) 
+ Krd{SlTi){SiT2) + K^{S2T2){t^S2) 
+ {SiS2)[K^{SiTi) + K[{S2T2)\ 

+ {SiS2)[K2{riS2) + K'^{SiT2)] 

+ Ul{SlS2f{TlT2) + U2{SlTi){SlS2){S2T2) 

+ U^{t^S2){S^S2){S^T2)'Eo-1. (1) 

Here S and t denote spin-1 and spin-^ operators, re- 
spectively, and symmetric ordering of spin-1 operators is 
assumed wherever it is necessary. 

We construct different singlet MP wave functions in- 
terpolating between a few simple VBSpSlates, and use the 
technique of "optimal ground states" til to find a family 
of Hamiltonians for which those wave functions are exact 
ground states. Among the members of this family, the 
following interesting representatives are found: 

(i) "AKLT-like" models which differ from the AKLT 
model by a few additional terms; 

(ii) biquadratically coupled chains which do not con- 
tain bilinear exchange interactions between the 5 = 1 
and S — ^ legs; 

(iii) a model with purely bilinear interactions; 

(iv) multicritical models with infinitely degenerate sin- 
glet ground states (for finite systems the degeneracy is 
exponentially large). 

By the construction of our MP ansatz, all the ground 
states are dimerized, and since we choose the initial 
Hamiltonian to be translationally invariant, they are 
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spontaneously dimerized. This., is in line with the re- 
cent field-theoretical argumentEil that sufficiently strong 
biquadratic interaction can cause spontaneous dimeriza- 
tion in S — ladder, implying also a similar two-particle 
structure of the spectrum ( "absence of magnons" ) . [The 
elementary excitation of a spontaneously dimerized sys- 
tem is a pair of domain walls in the dimer order.] 

The paper is organized as follows: Sect. || explains 
the principles of constructing the matrix product ansatz 
and the procedure of finding the set of exact solutions, 
Sect. II] contains the most important results describing 



representative models as listed above, and Sect. IV gives 
a brief summary. For convenience, technical details and 
general solutions are listed in the Appendices, so that 
those readers who are not interested in the details may 



go directly to Sect. Ill 



here single arrows | t)i I i) indicate the spin-i states, 
and |-|-), |0), |— ) denote spin-1 states. 

The matrix elements of any irreducible tensor opera- 
tor (for definiteness, let us choose T^*) in a fixed basis, 
according to the Wigner-Eckart theorem, are given by 



T^'^{M',M) EE {JM\tI'^\J'M') 

= flj, {JM\kq,J'M') 



(5) 



where the reduced matrix element Tj j. 



does not de- 
pend on M, M', q and thus can be absorbed into the 
free parameters c'^'^'^'^ in (^). We need k to be half- 
integer, therefore necessarily J ^ J' and one arrives at a 
(2 J -f 1) X (2 J' + 1) non-square matrix. Further, we fix 
the choice of defining it as 



II. THE MATRIX PRODUCT ANSATZ AND 
OPTIMAL GROUND STATES CONSTRUCTION 



V = (-) 



k — q 



rpK 



(6) 



We construct the MP ansatz for the model from 
rectangular matrices and as follows: 



fo) = TMg^g^gf ■ ■ ■ 9^n^i92n) 



(2) 



where the matrix g^'^ contains spin states of the i-th 
rung only, and the total number of rungs is 2N. We de- 
mand that l^'o) is a global singlet, then according to the 
approach proposed in Ref. ^ this dictates the following 
structure of the elementary matrices g^'^: 



Am/ 



(3) 



" 2 ' 2 



V2 



1\) 



3 3 



L,R 



^3) 



Here T^''^ are in general arbitrary matrix representations 
of irreducible tensor operators T*^* transforming under 
rotations R according to the V*' representation of the 
rotation group: 

q' 

and lipXfi) are the rung wave functions with total spin A 
and its z-projection /x. The quantities a = c'-^^s) and 
b = c^2'2) are free parameters. We use the compact no- 
tation I -fr), I -D-) for the rung states having A = 4, and 
1 3), 1 3) for the states with spin A = |: 



2 ' 



|3> = 
ID = 



V3 



T). IS) 



OT), U) 



x/3 



|0i> 



-i), 



(4) 



For our problem one can choose, e.g., J = 1 and J' — ^, 
then we have 3x2 and 2x3 matrices for Tl and Tr, 
respectively, so that the dimension of matrix space coin- 
cides with the total number of states of one rung. 

More complicated basis for T^"^ can be chosen, e.g., 
one may "decompose" J = 1 into two Ji = J2 = | and 
define 

T^i'J^^M',MiM2) = (JiMiJ2M2|T^«|Ji2JAf) 

= ?J,Ji,J„J,. {Ji2JM\kq, JiM,, J2M2) , (7) 

where J12 denotes the eigenvalues of the operator (Ji + 
■ Then, combining Mi and M2 into one "composite" 
index, one gets 4x2 matrices for T^^ and 4x2 respec- 
tively for . Another difference with the previous case 
is that now for k — ^ there are two possibilities: J12 = 1 

or 0, and the set of matrices T^^ acquires an additional 
free parameter w: 



T, 



L.R 



T. 



L.R 



wT, 



ig,0 



L,R 



For the sake of brevity, we will further refer to those 
two construction as "2 x 3" and "2 x 4" MP ansatze; the 
explicit form of the matrices we used can be found in the 
Appendix |A| [see Eqs. (Al) and ([A2[)]. One may think of 
the above construction as of interpolating between sev- 
eral VBS states shown in Fig. |^. It is easy to check, for 
example, that the 2x4 ansatz at {b/a = \/2, w = 0) 
and 2x3 ansatz at b/a = \f2 correspond to the same 
state with completely dimerized legs shown in Fig. |2a, 
another choice of parameters (ajb = —2^2, w ±a73) 
leads to the state with S —\ leg dimerized and 5=1 leg 
in the AKLT-type VBS state (Fig. |b), and, finally, the 
combinations [w = ±1/a/3, 6 = 0) and (6 = 0, a ^ 0, 
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w oc 1/a) give two "U-shape" dimerized states shown in 
Fig. |2p,d respectively. 

Our procedure of constructing the exact ground states 
foUows the ideas presented in Refs. |l^,|l^,^: we require 
the MP wave function (H) to be a zero-energy ground- 
state of the local Hamiltonian hi^i^i, which ensures that 
it is a ground state of the global Hamiltonian H (an op- 
timal ground state, in the terminology of Ref. |2l| ). This 
yields the following conditions: 

(i) hi^i+i annihilates all states being matrix elements 
of the two products g[gf-^i and gf gti-i- 



If one or more of Aj"' is zero, this may indicate an addi- 
tional degeneracy of the ground state. 

Eqs. #, (I), (i, (0) and ®, ®, (0) wiU be the 
basis for the further analysis. Now we proceed to con- 
sidering specific solutions of those equations in various 
particular cases. 



hMg^)=0, h,2ig^g^)=0; 



(8) 



(ii) all other eigenstates of hi2 have the energy e > 0. 
Then |^'o) is the zero-energy ground state of H; if one 
drops the constant term —Eq-1 in (|l|), the remaining 
Hamiltonian has the energy density Eq per rung. 

For further^ treatment it is convenient to write the 
Hamiltonian h in terms of projectors on the states l^^jl^) 
of the two-rung plaquette + with fixed angular mo- 
mentum. The complete set of the plaquette states [see 
Eqs. ([A^)] contains one multiplet with J = 3, three 
quintuplets, four triplets and two singlets, and thus one 
obtains: 



III. MODELS WITH EXACT GROUND STATES 

In this section we present a number of models with ex- 
act ground states being the most simple representatives 
of different classes of solutions mentioned in the Intro- 
duction. 



A. 2 X 3 MP ansatz 



1. General case 



k.l=l,2 



(fe,0|,T,(fc)\ 



(0 I 

Ml 



k.l=l.A 



A3E|*3m)(* 



3M| 



(9) 



fc,i=1..3 



M 



M 



Here obviously Aj^''^ 



(Aj'*^^)* because of the hermi- 



tian property of h. The complete set of the plaquette 
states can be divided into two subsets: local ground states 
which are contained in the matrix 



9,{k)\ I. _ 



JM 



, fc = l. 



I* 

products g^g^ and g^g^, and local eigenstates 



JM I 



which do not enter there. 



The conditions (i) mean that the local Hamiltonian h 
should project only onto the states l^'j^/'), and the mul- 



jIj^) have to be absent in Eq. 
in the following system of equations: 



tiplets \^ JM 



This results 



It is easy to verify that for the 2x3 MP ansatz the two 
matrix products g^g^ and g^g^ contain generally only 
the following multiplets: 



1*^ 
I* 



as ^ - i,T,(2) \ 



2M/ ' 



then the remaining multiplets can be chosen as 



I* 



IM ~ V'im)! I*2m) 



(12) 



,13 N 
2M) j 



(13) 



A 



A' 



(s,(fc);9.(0) 
J 

(e,(fc');s,(/c)) 



0, fc=l. 
= 0, fc= 1. 



fc' = 1, 



Ag) 



(10) 



which is essentially a system of linear equations in the 
Hamiltonian coupling constants J.., K , U,, [see Eq. (^]. 
The conditions (ii) require that all the eigenstates 



within the subspace determined by the basis I^'/a?) 
positive energy, which yields the inequalities 

Xf^>0, a=l... , 

where Aj"^ denotes the eigenvalues of the matrix 



have 



(11) 



The conditions (^lCl,llh now take the form 



>(2,2) 
^0 



= A, 
A 

= A 



(1,2) 


(fe,o 
1 

(2,3) 



= K 

= 0, 

= A 



(2,1) 



= 0, 

2<fc<4, A:<;<4, 



(3,3) 



= 0, 



(14) 



A^'"^ > 0, X\' '' > 0, A3 > 0, A^ > 0, 



where Aj , 01 — 1, 2 are the eigenvalues of the matrix 



5,(1,1) ,(1,2) 
^2 '^2 



1,(2,1) ,(2,2) 
■^2 ^2 
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(1 2) 

For the sake of simplicity we now set = 0. In 

fact, one can show that this requirement just fixes cer- 
tain "natural" symmetries in (hi), namely, 



Jr = J'r, Jd — J'd, Krd — K'^^, Ki,2 = K[ 



(15) 



Further, we require the six-spin couplings in (|^) to be 
zero, in order to make the model less cumbersome. Then 
one obtains the sytem of nineteen equations [fifteen equa- 
tions (|4|) and four additional assumptions X':^'^^ = 0, 
Ui.2,3 — 0] for twenty parameters of the Hamiltonian (|l|) 
and the free parameter u = a/b entering the 2x3 MP 
ansatz (|), (|). It turns out that one of those nineteen 
equations is linearly dependent, and the general solution 
contains, beside u, two additional free parameters x and 
y. This solution in its general form is presented in Ap- 
pendix ^ [see Eqs. (pl|), (B2)], and here we will just 
consider its most interesting particular cases: 

(a) Setting u — x = ||y and fixing the en- 

ergy scale by choosing 2/ = |, one gets the "AKLT-type" 
model of the form 



H — SnSn+l + - {SnSn+lY 



(16) 



+ 2 (S'nS'n+i) [(5,1 + Sn+l) ' (t n + T n+l)] , 

with the energy density per rung Eq = —2/3. (We recall 
that symmetric ordering of spin-1 operators is implicitly 
assumed). For this model two more eigenvalues of the 
local Hamiltonian are zero: 

Ao = Ai = , 

which may in principle indicate higher degeneracy of the 
ground state (there always exists another dimerized sin- 
glet state which can be obtained by the translation in 
one rung, and Ai = may mean degeneracy with some 
"partially ferromagnetic" state with the total spin J = 1 
of each plaquette) . 

(b) Setting u = —^^/2, y = 0, and choosing x — ^ to 
fix the scale, one obtains a model of S" = 1 and S ~ ^ 
chains coupled with purely biquadratic interaction: 



(17) 



with the energy per rung Eq = —2. For this model also 
one of the local Hamiltonian eigenvalues vanishes, 

A^^)=0. 

Another model of this type can be obtained by setting 
u = y = 0; after choosing a; = ^ to fix a proper 

scale the Hamiltonian takes the form 



+ -{SnSn+l) [{Sn + Sn+l) ' (r„ + T„+i)] (18) 
— [{Sn + Sn+l) ■ T n\ [{S n + Sn+l) ' T ri+l] ■ 

The energy per rung is i?o = —1, and the eigenvalue 
Ao = 0. 

(c) If a; = 0, then one can somewhat surprisingly ob- 
serve that the solution (Bl) does not depend on the 
parameter u = a/b entering the MP wavefunction (|^). 
Three of the local Hamiltonian eigenvalues are now zero. 



Ao = 0, A 



(1,1) 



0, A, 



0, 



and one has to put y < for the remaining ones to be 
positive. Setting the energy scale by fixing y — one 
obtains the model 



H — r„T n+l + -^SnSn+1 



(19) 



-S'„(2t n + T n-1 + T n+l) + Ij^iS nS n+l) {t nT n+l) 



1 



[{Sn + Sn+l) ■ T n] [{Sn + Sn+l) ' T n+l] 



with the energy density Eq = —j per rung. The fact 
that A3 = means that the singlet ground state of the 
model is degenerate with the fully polarized ferromag- 
netic state: the ferromagnetic state is the eigenstate of 
the Hamiltonian, and it is a straightforward exercise to 
check that it has the same energy. However, the ground 
state degeneracy of the model (|l9|) is much higher: any 
wave function \^{u)) of the form (||), (||) with arbitrary 
parameter u is the ground state. One can easily calcu- 
late the overlap between two such wave functions having 
different values of u: 



{^{u,)\^{u2)) 



N 



(l + muz)^ 

{l + ul){l + ul) 



<1, (20) 



i.e., the two g.s. wave functions with different values of 
u are asymptotically orthogonal in thermodynamic limit 
N 00 with the overlap vanishing exponentially with 
the increase of N. This means that the dimension of the 
basis of this subspace {[^'(w))}, i.e., the number of mu- 
tually orthogonal ground states, is exponentially large in 
thermodynamic limit. This is an example of multicriti- 
cal model. Unfortunately, within the present approach 
one cannot make any statement about presence of the 
gap. It is worthwhile to mention that infinitely degen- 
erate ground state in a mixed spin-1 /spin- ^ system was 
observed in certain limiting case of de Vega-Woynarovich 
model (see the discussion of c = case in Ref. I) when 
the velocity of one of the two spinon branches with linear 
dispersion becomes zero; however, since a pair of spinous 
can be combined either in a singlet or in a triplet with 
the same energy, the set of degenerate ground states in 
that case should contain not only singlets, but states of 
higher spin as well. 
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2. Case & = (no spin-^ states on the rung) 

For 2x3 MP ansatz it is useful to consider separately 
the case 6 = 0, because it turns out to lead to a new type 
of solution. One can see that at 6 = the spin-| states 
of the ladder rungs are excluded from the wavefunction 
and only the following two multiplets are present in 
the matrix products gig^, 9i'92- 



Mm) 



IMI 



(21) 



It turns out that in this case it is possible to obtain a 
nontrivial solution of the system ( p^ ) with only bilinear 
coupling, which correspond to the following model: 

H — ^ jSnT n — T nT n+1 — SnSn+1 

n 

- S„(r„_i +T„+i). (22) 

Here 7 > | is an arbitrary parameter, and the energy 
density is i?o = — (7 + |) per rung. At 7 = | the 
eigenvalue A3 vanishes, indicating the first-order transi- 
tion into fully polarized ferromagnetic state [cf. a similar 
transition for S — ^ ladder in Rcf . ^ . The matrix prod- 
uct gig2 has the simple form 



2| u) - 1 m 



U) + 1 U) 



(23) 



which allows one to "visualize" the structure of the 
ground state as interpolating between two "U-shape" 
VBS states shown in Fig. |;,d [here i[, 4 are the rung 
states with total spin i, see Eqs. (^)]. 



them, yielding a one-parametric family of multicritical 
models with infinitely degenerate ground state similar to 
one discussed in the previous subsection. The solution is 
given by 

J, = j; = 2x, Jd = J'd = Ki=K[=K2 = K'2^x, 
3 4 

Ks^\v. i^sr = ^(132/- 9x), (26) 

4 

Krr = Kdd = 3a; - 2y, Krd = K'^d = y(3a; - 2y) , 
;7i = ^(5y-4a;), U2 ^ = ^{ix - 2y) , 

and the conditions ( |ll| ) take the form 

20 

A2 = 8(y - .x) > 0, A3 = y(2x + y) >0. (27) 

This class of solutions is more cumbersome than in case 
of the 2x3 MP ansatz; for instance, one may observe 
that six-spin interactions C/1,2,3 in ( |2^ ) should be always 
nonzero. The simplest model within this class is achieved 
by setting y = |x, its Hamiltonian after fixing the proper 
energy scale is 

H — ^ SnSn+l + t: [SnSn+lf [l + TtT nT n+l] 



gT nT n+l + g S'„ (2t n + T„_i + T n+l) (28) 
-{SnSn+l) \{Sn + Sn+l) ' (t„ + T„+i)] , 



B. 2 X 4 MP ansatz 



For the 2x4 MP ansatz, at a general choice of param- 
eters a, &, w entering the wave function (|^), the following 
multiplets are contained in the two matrix products g^g^ 
and g^g^: 



with the energy per rung Eq = It should be re- 

marked that, in contrast to the model (|l9|), fully polar- 
ized ferromagnetic state generally (except for the case 
X — —■^y) is not degenerate with the ground state of 
multicritical models defined by (p6|). 



(1/V2)(V. 



2M 



2m)^ '^2M 



13 



(24) 



We require them to be annihilated by the local Hamilto- 
nian h, and the remaining multiplets 



V3M 



(25) 



to be the eigenstates of h with positive energy. Eqs. (|1C 
give a system of eighteen linear equations for twenty pa- 
rameters of the Hamiltonian (||), and its general solution 
contains two free parameters x, y (one of them is again 
irrelevant since it just sets the energy scale). The multi- 
plets (|2^) do not contain the wave function parameters 
a, 6, w, and thus the solution also does not depend on 



IV. SUMMARY 

We have studied the isotropic ladder composed of 
5 = 1 and S = ^ chains, with general type exchange in- 
teraction between spins on neighboring rungs. The tech- 
nique of matrix product states is applied to construct a 
family of Hamiltonians with exact ground states. Among 
the members of this family, we have found a couple of 
relatively simple models with nontrivial ground states, 
including one model with only bilinear exchange and two 
models with S = 1 and S = ^ chains being coupled by 
purely biquadratic exchange. We also present a family 
of multicritical models whose ground state is infinitely 
degenerate: in the thermodynamic limit N —^ oo the 
number of degenerate ground states grows with N expo- 
nentially. 
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APPENDIX A: TECHNICAL DETAILS 

Here is the explicit form of the matrices we used in (|^) 
for the 2 X 3 MP ansatz: 



^-3^3 = (1/2V5){|33) + |33) + 3(|11) + |11))}, 



^ L 



3 3 
7^ 2 ■ 2 

— 



7^ 2 ■ 2 







1 
















1 


2 




rri 2 ' 


3 _ 








3 


3 


Tl 


2 _ 











T 


fi 






V 3 

















1 



(Al) 















The matrices Tfj can be easily obtained from those ma- 
trices using the definition (^. 

And in case of the 2x4 ansatz the matrices were chosen 



as 



1 


1,0^ 


1 


-1,0 


1 




1 

2 




3 

7^2 


3 

2 


3 

7^2 


1 

2 


3 

7^ 2 


1 

2 



(A2) 



1/V2 -l/\/2 




1/V2 -l/\/2 

1/V6 l/\/6 
^2/3 

-72/3 

-l/\/6 -l/%/6 

l" 


-i/Vs -i/Vs 

1/V3 

1/V3 

-l/\/3 -l/\/3 



The complete set of spin states of a two-rung plaquette 
is, in the notation used in Eqs. (0), 



7^2 ' 2 
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V'33-|33>, ^'3''-3 - |33), 

^i = i=(|31)-f |13)), V'|'-2 = ^ 

V'l'i-(l/V5)(V3|ll) + |31) + |i3)), 
V'!'-i = (l/V5)(\^|ll>-K|31) + |13)), 



(|31) + |13)), 



Vi = -^(|31>-|13)), 



1 

71^ 



= ^(|3i) - |13)), ^f_. 



1 

1 

71 



(|3T) - |13)), 
(|31)-|13)), 



^2^i^ = (l/2)(|ll)-|ll) + |33)-|33)), 

= VWlO(|31) + |T3)-(2/V3)|ll)), 
V'f.i = 73710(131) + |13) - (2/^/3)|lT)), 

= (l/2y5){3(|33) + 133)) - |lT) - |11)}, (A3) 
V^o'3 = (l/2)(|ll)-|ll) + |33)-|33)), 

^2^2 = 13^), ^2'-2-|3^>, 

i^ll = ^|3 ^) + ^|1 ^), V2' -1 = ^|3 ^> + 4), 



1 



V'2'0 = ^) + |1 - ^> + ^|3 ^), 



V3, 
2 



^f-i - olU) - ^|3 ^), i^ll = ^(|i 4) - |i ^)), 



^1 

1 

Vl^ 
1 



V3, 
2 



C = ^(1 n) - 1 M)), i^ll - ^(1 M) + 1 M)), 

^n = IM>,^i\-i-|M). 

Here the superscripts of V''s denote the total momentum 
J of left and right rung states of which they are com- 
posed (with the shortcut convention that 3 means J — ^ 
and 1 means \). The states V'jm can be obtained from 
V'^m interchanging the left and right rung states. 



APPENDIX B: GENERAL SOLUTION FOR THE 
2x3 MP ANSATZ 

The general solution of the system dl^) with four addi- 

(1 2) — 
tional conditions A2 — 0, ?7i_2,3 = has the following 

form: 

Js = [iSu" -f 371*2 ^y(2u4 - 40^2 - 3)] a; - 2/ , 

Jr = (112U-* + 176u^\/2 + 192^2 + uV2u + U)x - 8y , 

Jr = 2V2(4u + V2){u^ + l){2u^ - l)x + 2y , 

Jd^y, (Bla) 

Ks ~-V2{2u + ^/2){u + V2)u(-2u -I- V2fx , 
Ksr =4[13?/ + 2511^ + 3/2 

- V2m(2u'* - 24^2 - 5)] X - 4y , 
Krr = -V2{-4:U* + llu^V2 + 2\v? + 7\/2m -|- 5) 

X (2u + \/2)a; + 4y , 
Kdd = - V2(4w'' 9m^V2 + 271*2 -I- 5\/2u + 3) 

X (2u -f ^Pi)x -f 4y , 
-ft^rd = -4(5u2 + \/2u + l)(w + \/2) 
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X (2u + \/2)x + 4y , 
Ki = -{-2u + V2){u + \/2)(u2 + 3\/2?i + l)x , 



V2) {8u^ + 5m^\/2 + 2u - V2)ux , 



K2 = i\/2(-2M 
Ea = -[44^^^ + 68^2 + 9/2 

+\/2u(8u'* + eOu^ + 13)] X + 2y ; 

with the conditions on the eigenvalues ( p+j ) being 
Ao = -6\/2(4?i + ^/2)x > , 

Ai = (48u'* + 66u^\/2 + 78^2 + 18\/2u + 6)a; - 4y > , 
A^^'^^ =3\/2(2M + \/2)(u2^1)a;>0, (Bib) 
Xl'^ = 12(121*3 ^ 9^272 + 4u + \/2) 

x(u + \/2)x-12y>0, 

A3 =60u2(m + \/2)2.t > 0. 

Here u = a/6 is the free parameter entering the 2x3 
MP ansatz (^, and x,y are additional free parameters 
arising from the solution of the linear system ([T^). Tak- 
ing into account that one of the parameters just se ts t he 
energy scale and thus is irrelevant, one gets from (Bl) a 
two-parameter family of Hamiltonians with exact ground 
states in a form of the 2x3 matrix product. One can see 
that inequalities ( Blb| ) can be satisfied provided that 



\/2 V2 

x> 0, <u< - — 

- 4 - - 2 



(B2) 
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FIG. 1. Mixed spin ladder with the S = i upper leg and S = 1 lower leg as described by the model (jl]) 




(a) 




(b) 




(c) 



id) 



FIG. 2. A couple of VBS states which can be obtained as limiting cases of the 2x4 and 2x3 MP ansatze (^). Dotted box 
corresponds to the matrix product g^g^. 
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